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As one of the key units of space CCD camera, the temperature range and stability of CCD components affect the 
image’s indexes. Reasonable thermal design and robust thermal control devices are needed. One kind of tempera- 
ture control loop heat pipe (TCLHP) is designed, which highly meets the thermal control requirements of CCD 
components. In order to study the dynamic behaviors of heat and mass transfer of TCLHP, particularly in the or- 
bital flight case, a transient numerical model is developed by using the well-established empirical correlations for 
flow models within three dimensional thermal modeling. The temperature control principle and details of mathe- 
matical model are presented. The model is used to study operating state, flow and heat characteristics based upon 
the analyses of variations of temperature, pressure and quality under different operating modes and external heat 
flux variations. The results indicate that TCLHP can satisfy the thermal control requirements of CCD components 
well, and always ensure good temperature stability and uniformity. By comparison between flight data and simu- 
lated results, it is found that the model is to be accurate to within 1°C. The model can be better used for predicting 


and understanding the transient performance of TCLHP. 
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Space CCD camera 


Introduction 


Space CCD camera has been used in a wide range of 
fields, including military reconnaissance, resources ex- 
ploration, surveying and mapping etc. [1] As one of the 
key units, the temperature range and stability of CCD 
components affect the image’s indexes of signal-to-noise 
ratio and geometric accuracy directly or indirectly [2-4]. 
For some high performance camera, traditional products, 
limited by the heat transfer capacity, mounting dimen- 
sions and temperature stability, are not suitable for the 
thermal control of CCD components [5]. Reasonable 
thermal design and robust thermal control devices are 
needed. 


Nomenclature 
Symbols 
A area (m^) 


c specific heat at constant pressure (J/(kg:°C)) 


p 


LHP is an advanced heat transfer device, which uses 
the capillary pressure generated in a porous structure to 
circulate the working fluid from heat source to heat sink 
[6-7]. For the purpose of precise thermal control of CCD 
components, one kind of TCLHP is designed. By de- 
coupling the evaporator and heat source, the capillary 
force is generated by loading heat directly on the evapo- 
rator, which force drives the fluid to absorb the heat pro- 
duced by CCD components and release it through con- 
densers. Compared with the traditional use of LHP, 
TCLHP has the characteristics of well start-up, good sta- 
bility and flexible piping layout, which is promising in 
the thermal control of CCD components. 

Extensive ground and flight tests have demonstrated 


Greek symbols 
A difference 
y latent heat of vaporization (kJ/kg) 


D diameter (m) u dynamic viscosity (Pa's) 
G heat conductivity (W/K) p density (kg/m°) 
h convective heat transfer coefficient (W/(m?:°C)) Subscripts/superscripts 

k heat conductivity (W/(m-°C)) accu accumulator 

K permeability (m°) cw capillary wick 
L length (m) e evaporator 

m mass flow (kg/s) env environment 
Nu Nusselt number f fluid 

P Pressure (N/m°) i inner 

Pr Prandtl number l liquid 

Q heat load (W) o outer 

r radius (m) sat saturated 

Re Reynolds number vg vapor groove 
T temperature (°C) w wall 

u velocity (m/s) 

x quality 


the excellent heat transfer performance of LHP [8]. 
However, compared with the experimental investigations, 
only a few studies on the theoretical analyses and ma- 
thematical modeling of LHP have been conducted so far. 
Accurate LHP model is needed to study the thermal per- 
formance, particularly in the orbital transients. The LHP 
system involves complex heat and mass transfer pro- 
cesses, rendering the mathematical modeling highly 
challenging. A number of steady-state models with vary- 
ing levels of complexity have been proposed, which 
prove to be useful in predicating steady-state perfor- 
mance characteristics. Kaya et al. [9] developed a one 
dimensional (1D) mathematical model of LHP, which 
could reflect the variable conductance characteristics, but 
the oversimplification in the model brought large differ- 
ence between the modeling and experimental results. A 
1D steady-state model based on the energy balance equa- 
tion, thermodynamic relationships and heat transfer and 
pressure drop calculations was developed by Chuang 
[10]. Improvement included the pressure drop induced by 
the bends, the convective heat transfer in the vapor 
grooves and heat fluxes in the wick. Adoni et al. [11] 
established a model to study the thermal and hydraulic 
performance of LHP, by using the mass and energy equa- 
tions in the system. Additionally, the pressure drop across 
the capillary wick is calculated by the Darcy’s law. Bai et 
al. [12] developed a model to study the effect of a 
two-layer compound wick, based on the energy conser- 
vation equations and empirical correlations. Rivière et al. 
[13] modeled a LHP based on the classic nodal network 
to study the fluid distribution in a loop. 

Only a few of mathematical models have been estab- 
lished to study the transient behaviors in response to the 


power or sink cycles. For example, Wrenn et al. [14] and 
Hoang and Ku [15] developed a transient model based on 
1D time-dependent conservation equations, which was an 
extension of transient model of capillary pumped loop 
(CPL). There were large errors between simulation and 
experimental results. Pouzet et al. [16] presented a com- 
prehensive mathematical model, where the whole dy- 
namic behaviors of a CPL due to power steps were suc- 
cessfully captured. Their model is helpful for the tran- 
sient simulation of a LHP. Cullimore and Baumann [17] 
developed an advanced node-type LHP model written by 
FORTRAN code, which was embedded into Sinda/Fluint 
analyzer. This model is used by several researchers, such 
as in Refs. [18] and [19]. Compared with the test results, 
the predicted temperature showed little differences (less 
than 1°C). For the aspect of orbital simulation, Ferrandi 
et al. [20] established a transient model to design and 
simulate the thermal control system of LHP by Sin- 
da/Fluint analyzer, which is used for the cooling of the 
cryomagnet avionics box on the Alpha Magnetic Spec- 
trometer (AMS-02). Xin et al. [21] built a system LHP 
level for the thermal control of cryocooler on AMS-02, to 
analyze the transient performances due to the orbital 
change environment, and to understand the operation of 
LHP during the mission time. In any case, it is capable of 
aiding system design and predicating the working per- 
formance on orbit by LHP models based on the Sin- 
da/Fluint code. 

In this study, a transient mathematical model of 
TCLHP is established by using Sinda/Fluint code. In the 
following, the operating principle and mathematical 
model are presented first, followed by numerical results 
and discussions and comparison between flight data and 


simulated results. Finally, some conclusions are drawn. 
Analysis and Modeling 


Operating principle 

For the purpose of heat dissipation and temperature 
control of multi heat sources on the space CCD camera, 
the flow schematic of TCLHP is shown in Fig. 1. The 
operating principle is as follows. As heat is applied to the 
evaporator, liquid is vaporized and the capillary force is 
developed to push the vapor out of the evaporator. Vapor 
condenses into liquid in the condenser A and the capillary 
force continues to push liquid into the pre-heater, in 
which the fluid is changed into two phase fluid. Then the 
saturated fluid is forced into the cold plates and heat 
produced by CCD components is absorbed by the latent 
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Fig.1 The flow schematic of TCLHP 


heat while the temperature keeps near constant. Two- 
phase fluid condenses into sub-cooled liquid again in the 
condenser B and the capillary force pushes it into the 
accumulator. And that cycle repeats. 

Because the saturation states of fluid that in cold 
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plates and accumulator are related, the following condi- 
tion must be satisfied [8]: 


AP, —accu 

CCD (1) 
(dP/dT ) 
where Tccp and Taccu are the saturation temperature of the 
fluid in cold plates and accumulator, APocp-accu 18 the 


Tecp ~ Laces = 
sat 


saturation pressure difference between cold plate and 
accumulator, and (dP/dT),q, is the slope of the pres- 
sure-temperature saturation line at T,,.,. Eq. (1) states 
that, for a given pressure difference between these two 
elements, a saturation temperature difference corres- 
ponding must also exist. Generally, the number values of 
(dP/dT),q are relative large. For example, the value of 
(dP/dT)sa at 5°C is 3.78 kPa/°C. For TCLHP, the mass 
flux is small and the tubes are smooth, so the pressure 
difference is small. If the number value of APocp-accu 


is 1 kPa, the temperature difference is only 0.26°C ac- 
cording to Eq. (1). Hence, the working temperature of 
CCD components can be controlled by stabilizing the 
temperature of accumulator. 


Mathematical Model 


The evaporator and accumulator 


The schematics and thermal network of evaporator 
and accumulator are shown in Fig. 2. Since the accumu- 
lator is integrated with the evaporator, there are compli- 
cated relations of heat and mass transfer between them. 

The heat loaded on the evaporator wall is divided into 
three parts. The first part is from evaporator wall to the 
wet wick in evaporator by convective heat transfer, as 
shown by 
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Fig. 2 The schematics of evaporator and accumulator (a), and the thermal network (b) 


Om =R-(th,DeL,)-(Tye-Ty¢) (2) 
with R an empirical constant that accounts for the heat 


transfer in the wick. The resulting value for R is a com- 
bination of conduction, convection, boiling, and evapora- 
tion between the wet wick and the evaporator wall [22]. 


This part also splits into an evaporation power and the 
back-conduction into the accumulator. 
The evaporation power can be calculated by 
Qn =Y (3) 
The back-conduction power from wet wick to accu- 
mulator is given by [12] 
mc 


Oini2 = (D a y a x (Ty. T accu) (4) 
where 
MC p 
1” rk La ©) 


This section of part is related with the properties of 
capillary wick and working fluid and the mass flux. 

The other two parts of heat are radial heat leak from 
the evaporator to the accumulator and the heat used for 
the superheating the vapor in grooves. Thus, 


Qina T Suas (Toe ~ Ty ,acou ) (6) 
Oin3 T h, A (Tre = Tove) (7) 


The thermal balance equation of accumulator satisfies 
the following relation: 


Oe oF Onz + 0.9 z Osy + Qub (8) 
where Qaccu is the power of temperature control loaded 
on the accumulator, and Qeny and Qsup can be calculated 
by 

Oiv = Gen (Tp accu 7 To) (9) 


Qsub = MC py ea, z Ty sub) (10) 


Transport tubes 

The transport tubes of TCLHP are usually small di- 
ameter smooth pipes with long transport length. By neg- 
lecting the heat exchange between tubes and environment, 
the energy conservation equation of the working fluid in 
lines can be written as [12]: 


=j 
. dT 1 In(D 1D:) 
= th 11 
Cpl dL i 2mk ( J onv) ( ) 
The heat transfer coefficient at the inner wall can be 
calculated by 


cw 


Nu )k 
h, = Wu)ks (12) 
D, 
Since there is large length-to-diameter ratio for the 
transport lines, the developed flow is assumed. Nu can be 


determined by the following equations: 


4.36 Re < 1960 
(f /8)(Re-1000) Pr 


2 1960 < Re < 6420 (13) 
1+12.7(f/8) (Pr?-1) 


Nu= 


0.023 Re®® pr” Re > 6420 


where n is 0.4 when the pipe wall heats the fluid and n is 
0.3 when the pipe wall cools the fluid. 

The pressure drop in the transport lines can be calcu- 
lated by 


dP l 3 4 
L py x 14 
A (14) 
where the frictional factor is calculated as follows: 
64/Re Re < 2200 (15) 
~ |0.0791Re°?>  2200<Re < 10° 


Two-phase flow model 


Two-phase flow phenomena exist in pre-heater, CCD 
cold plate and condensers. The heat transfer coefficient 
of two-phase is related with the heat flux density, stabili- 
ty and flow pattern. 

Shah correlation is used to calculate Nu as follows 
[23]: 


0.8 p,.0.4 og 3.8x°" (1 -x)™ 
Nu =0.023 Ref’ Pr?4| (1—x) "+ SOR (16) 
where 
Nu= ve (17) 
Hc 
Pr, = ie (18) 
P= 2 (19) 


The qualitative temperature in Eqs. (16)-(19) adopts 
the saturation temperature of two-phase section. 


Capillary wick and vapor grooves 

According to Darcy’s law [24], the pressure drop 

through the capillary wick is determined by 
oar (20) 
271K ewi Lewi 

The vapor is generated along the whole groove length, 
its mass flowrate increases approximately linearly along 
the vapor grooves. Hence, the effective length of vapor 
groove is approximated as half of its total length. 

The pressure drop of vapor grooves can be calculated 
by Eqs. (14)-(15). In order to determine the Reynolds 
number, the hydraulic diameter of each vapor groove is 
used: 

2tw 


E T (21) 


where ¢ and w are the height and width of the vapor 
groove, respectively. 


Parameters and Conditions 


Model Parameters 


The working fluid is ammonia. The temperature of 
accumulator wall is controlled at 5°C during simulations. 
Table 1 shows the main parameters of each component. 


The external heat flux variations and the working 
mode of CCD components 
Fig. 3(a) shows the simulated external heat flux varia- 


tions of condenser A and B. Fig. 3(b) gives the working 
mode of CCD components. The orbital period is 5852 s, 


Table 1 Model parameters 


Components Description 
. Material: ceramic; Pore radius: 1x10% m; 
Capillary wick : 
Porosity: 0.6 
Material: stainless steel; Length: 0.13 m; 
Evaporator 
OD: 0.014 m; Channels: 20 
Material: stainless steel; 
Accumulator i 
Volume: 1.1107 m 
Material: aluminum alloy; Area: 0.35 m° for 
Condenser 3 
condenser A, 0.33 m” for condenser B 
Material: stainless steel; 
Tubes 


OD: 0.003 m; ID: 0.002 m 
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Fig.3 The external heat flux variations of condenser A and B 
(a); the working mode of CCD components (b) 
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and the operation period is 900 s. During operation, the 
power dissipation of each CCD component is 9 W, and 
the total heat power is 36 W. 


Results and Discussion 


The temperature contour maps of condenser A and B, 
for cold case (a) and (b), and for hot case (c) and (d), are 
shown in Fig. 4. The cold case is corresponding to the 
time instant of the smallest external heat flux at the cam- 
era shut-down, and the hot case is corresponding to the 
time instant of the largest external heat flux during oper- 
ation. In the application, the vapor line is coiled on the 
upper right section of the condenser B firstly, and then 
enters into the condenser A at which the fluid is sub- 
cooled enough. The working fluid is heated to be satu- 
rated in pre-heater. The two-phase fluid absorbs the heat 
produced by CCD components. The “hot” pipe is then 
coiled on the lower left of condenser B. After the work- 
ing fluid changes to be sub-cooled again, it flows into the 
accumulator. The bulk temperature of condenser A and B 
is low in the cold case, while it is high in the hot case. 
With comparison of the four contour maps, it can be ob- 
served that the sections on which vapor and two-phase 
lines coiled are hot. 

Fig. 5(a) and (b) show the temporal evolution of tem- 
perature and quality of working fluid in the pre-heater. 
The pre-heater is used to heat the sub-cooled liquid to 
saturated condition (i.e. set-point) before entering into 
the CCD cold plates. The numbers from 1 to 15 represent 
the fluid at different locations from the inlet to the outlet 
of pre-heater. Fig. 5(a) indicates that the external heat 
flux and working mode of CCD components affect the 
temperature of fluid in pre-heater. The temperature of 
outlet fluid is near 4.5°C after heated in pre-heater. Fig. 
5(b) shows that the quality changes from 0 to 0.12-0.17, 
which means that the state of working fluid has changed 
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Fig. 4 The temperature contour maps of condenser A and B, 
for cold case (a) and (b), and for hot case (c) and (d) 
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Fig. 5 The temporal evolution of temperature (a) and quality 
(b) of working fluid in pre-heater 


to be saturated. The pre-heater assures that the fluid 
flowing into CCD cold plate always keeps in two-phase 
condition during the mission time. 

Fig. 6(a) and (b) give the temporal evolution of tem- 
perature and quality of working fluid in CCD compo- 
nents. After heated in pre-heater, the fluid that entering 
the cold plate becomes two-phase saturation state. During 
the camera shut-down period, the saturated fluid flows 
though the CCD cold plates in a near constant tempera- 
ture. During operation period, the liquid in two-phase 
fluid is evaporated, and thus the quality increases. The 
temperature of fluid keeps in the range of 4.37-4.53°C, 
which provides a stabilized working boundary for CCD 
components. There is a decreasing temperature trend 
along the four cold plates due to the two-phase pressure 
drop. The temperature uniformity can be assured by the 
fluid characteristics of two-phase fluid. Unlike the tem- 
perature change, the quality changes obviously. The 
quality increases from 0.12 at the inlet of the first CCD 
cold plate to 0.90 at the outlet of the fourth CCD cold 
plate. The results indicate that the influences of external 
heat flux and working mode on temperature can be ad- 
justed by the quality of two-phase fluid. 

The profiles of temperature, pressure and quality of 
working fluid along flow distance, for cold case (a) and 
(b) and for hot case (c) and (d), are shown in Fig. 7. 
During the two cases, the temperature and quality de- 
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Fig. 6 The temporal evolution of temperature (a) and quality 
(b) of working fluid in CCD components 


crease firstly, corresponding to that the vapor from eva- 
porator becomes sub-cooled liquid. In order to save pow- 
er, heat exchange between the sub-cooled and return two- 
phase lines occurs. Then the liquid flows into the pre- 
heater, and it turns into two-phase fluid. The temperature 
increases to a saturated point and the quality increase 
from zero to a number greater than zero. In cold case, the 
temperature and quality of two-phase fluid keep un- 
changed along the four CCD cold plates. In hot case, the 
temperature stays the same while the quality increases 
quickly along the heating elements. Saturated two-phase 
fluid condenses into sub-cooled liquid, and then it is 
pushed into the accumulator. The pressure is dropping 
along the flow direction. 

Fig. 8 shows the comparison of average temperature 
of four CCD cold plates for the states of non-operation 
period (a) and operation period (b), between the data in 
orbit and the simulation results. During shut down period, 
CCD cold plates keep near 4.7°C in orbit, and the simu- 
lated temperature is around 4.4°C. The temperature dif- 
ference between tests and simulation results is only 0.3°C. 
During operation, the average temperature of four CCD 
cold plates increases from 4.9°C to 5.5°C for flight re- 
sults, and the temperature of them increases from 4.8°C 
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Fig. 7 Profiles of temperature, pressure and quality of working fluid along flow distance, for cold case (a) and (b) and for hot case 
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Fig. 8 The comparison of average temperature of four CCD cold plates for the states of non-operation period (a) and operation pe- 


riod (b) 


to 4.9°C for simulation results. The temperature differ- 
ence between the two cases is in the range of 0.1°C ~- 
0.7°C. The model slightly under-predicts the loop oper- 
ating temperature. The error between the model and the 
curve fit is generally less than 1°C, in many cases, less 
than 0.5°C, which verifies the validity and reasonability 
of TCLHP model. The model can be used for predicting 
the transient characteristics of LHP, and for understand- 
ing and analyzing the TCLHP’s performance. 


Conclusions 


Numerical simulations have been carried out for a 
better understanding of the phenomena of heat and mass 
transfer for TCLHP of space CCD camera. The main 
conclusions from the numerical simulations are given as 
follows: 

(1) The fluid flow states in the tubes coiled on the 
condensers affect the temperature distribution of con- 


denser A and B. The temperature of the sections of vapor 
and two-phase lines coiled on is high, while the other 
sections are cold. 

(2) The pre-heater between the condenser and CCD 
cold-plates assures that the fluid flowing into CCD cold- 
plate keeps in two-phase condition always during the 
mission time. 

(3) The temperature of fluid in cold plates is stable and 
uniform. The influence of external flux and working 
mode on temperature can be adjusted by the quality of 
two-phase fluid. 

(4) Comparison between flight data and simulated re- 
sults indicates that the model slightly under-predicts the 
loop operating temperature and the error between the 
model and the curve fit is generally less than 1°C, in 
many cases, less than 0.5°C. 

The results indicate that the temperature boundary of 
CCD components can be provided by the TCLHP, which 
can assure the CCD components working normally. The 
model can be better used for predicting and understand- 
ing the transient performance of TCLHP. 
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